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ABSTRACT 


At  electron  densities  greater  than  10^^/cc,  ionization  and  recombina¬ 
tion  proceed  predominantly  by  sequential  transitions  between  adjacent 
energy  levels  induced  by  collisions  with  free  electrons.  Similar  models 
for  this  process  have  been  developed  by  various  investigators.  This  report 
develops  closed-form  expressions  for  ionization  and  recombination  of 
effective  one -electron  atoms  by  modifying  the  Byron  model  and  determining 
analytically  the  critical  energy  level  in  equilibrium  with  the  free  electrons, 
using  a  simplified  quantum  mechanical  model  of  the  alkali  atoms.  The  re¬ 
action  rate  coefficients  calculated  from  these  expressions  agree  satisfac¬ 
torily  with  available  experimental  data  and  with  numerical  calculations  of 
other  investigators. 
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SECTION  I 
INTRODUCTION 


For  many  years  physicists  have  been  concerned  with  the  mecha¬ 
nisms  of  ionization  and  recombination  in  both  laboratory  and  stellar 
plasmas.  Recently  new  impetus  has  been  given  to  these  studies  by  two 
potential  applications  of  plasma  physics:  the  thermonuclear  program 
and  the  development  of  magnetohydrodynamic  (MHD)  energy  conversion 
techniques.  In  the  latter  case  research  has  concentrated  on  achieving 
conductivities  sufficient  for  practical  devices  by  seeding  plasmas  with 
small  amounts  of  alkali  metals.  As  a  result  electron  densities  of 
10^^  -  loi^/cc  can  be  attained  at  plasma  temperatures  low  enough  that 
conventional  gas  dynamic  channels  can  be  used  for  containment  of  the 
plasma.  Because  the  gas  velocity  in  such  channels  is  usually  quite  high 
{on  the  order  of  1000  m/sec),  there  are  conditions  under  which  a  sizable 
fraction  of  the  channel  length  is  required  for  a  perturbed  plasma  to  reach 
equilibrium  conditions.  Even  when  equilibrium  conditions  prevail  in  the 
center  of  the  channel,  one  normally  finds  that  equilibrium  conditions  are 
not  maintained  near  the  electrodes.  Under  these  conditions  it  is  neces¬ 
sary  to  know  the  ionization  and  recombination  rates  in  alkali- seeded 
plasmas  to  fully  understand  the  flow. 

There  are  various  mechanisms  of  recombination  in  plasmas: 
dissociative  recombination  of  molecular  ions,  radiative  recombination, 
dielectronic  recombination,  three-body  recombination  with  a  neutral 
particle  participating  as  the  third  body  to  conserve  momentum,  and  three- 
body  recombination  with  an  electron  acting  as  the  third  body.  The  asso¬ 
ciated  mechanisms  of  ionization  are  listed  respectively:  associative  atom- 
atom  collision-induced  ionization,  photoionization,  autoionization, 
atom- atom  or  atom-molecule  collision- induced  ionization,  and  electron- 
atom  or  electron-molecule  collision- induced  ionization. 

Although  different  mechanisms  prevail  in  atmospheric  plasmas, 
conditions  occurring  in  laboratory  plasmas  and  practical  MHD  devices 
are  dominated  by  collision-induced  ionization  and  recombination  with  proper 
account  taken  of  radiative  effects. 

Understanding  of  the  general  process  has  come  about  only  recently. 
Massey  and  Burhop  (Ref.  1)  indicate  that,  at  the  time  of  their  publication, 
recombination  in  high  density  afterglows  resembled,  but  exceeded  by 
several  orders  of  magnitude,  pure  radiative  recombination.  In  an  effort 
to  explain  the  anamolously  large  electron  loss  occurring  in  the  B-1 
Stellerator,  D'Angelo  (Ref.  2)  applied  the  theories  of  non-thermal  equilib¬ 
rium  in  stellar  atmospheres  developed  in  Refs.  3  and  4  to  the  calculation 
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of  three-body  (electron-electron-ion)  recombination  into  excited  states  of 
the  resulting  neutral  atoms.  He  assumed  that  these  excited  states  de¬ 
cayed  further  only  by  radiant  emission.  Using  the  available  transition 
probabilities,  D'Angelo  computed  the  total  rate  of  recombination  into  the 
ground  state  from  all  possible  combinations  of  three -body  recombination 
into  an  excited  state  followed  by  radiative  cascading  downward  from  that 
state.  His  results  were  much  larger  than  those  caused  by  radiative  re¬ 
combination  alone;  however,  he  underestimated  the  results  of  experiments 
reported  in  Refs.  5  and  6  by  about  an  order  of  magnitude  and.  more 
importantly,  obtained  an  incorrect  dependence  on  electron  temperature 
and  density. 

During  this  time  Bates,  Kingston,  and  McWhirter  in  several  publica¬ 
tions  developed  a  statistical  model  of  the  process,  culminating  in  their 
two  authoritative  papers  on  collisional-radiative  recombination  (Refs.  7 
and  8).  This  model,  which  gave  better  agreement  with  experiments  than 
D'Angelo's,  considers  the  interaction  of  radiative  recombination,  three- 
body  recombination  into  excited  states,  radiative  and  collisional  cascading 
between  states,  and  excitation  and  ionization  from  various  states.  Precise 
quantum  mechanical  expressions  were  used  for  the  various  radiative 
transistions,  and  the  semiclassical  cross  sections  of  Gryzinski  (Ref.  9) 
were  used  for  the  collisional  transitions.  Qualitative  information  on  the 
behavior  of  alkali  plasmas  was  obtained  by  making  the  firs't  excited  state 
of  hydrogen  the  ground  state  of  a  fictitious  substance  having  a  potential 
well  about  the  depth  of  that  of  cesium.  The  results  indicated  a  similarity 
of  behavior  of  hydrogen  and  alkali  plasmas. 

At  about  the  same  time,  Hinnov  and  Hirschberg  (Ref.  10)  developed  a 
theory  to  correlate  the  Stellerator  experiments  by  determining  separately 
collisional  and  radiative  recombination  coefficients.  In  their  paper 
collisional  transition  rates  between  all  states  up  to  some  state  in  equilib¬ 
rium  with  the  free  electrons  were  determined  by  using  Thomson-Bohr 
cross  sections.  The  lowermost  state  in  equilibrium  with  free  electrons, 
denoted  henceforth  as  the  critical- state,  was  located  by  assuming  detailed 
balance  for  excitation  and  de- excitation  of  this  state  by  collisions  and 
radiative  transitions.  The  resulting  relation  was  evaluated  in  the  limit  of 
high  critical  quantum  number.  For  electron  energies  of  0.  25  ev  or  less, 
the  approximate  collisional  recombination  coefficient  obtained  by  this 

/k  TV  ”4 

method  for  hydrogen  is  5.6  x  10”27  cc/sec.  Good  agreement 

with  experiments  in  both  hydrogen  and  helium  plasmas  was  claimed  when 
the  radiative  recombination  coefficient  was  added  to  the  results  quoted 
above. 
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Byron,  et  al.  (Ref.  11)  obtained  almost  identically  the  same  results 
for  hydrogen  at  high  temperature  as  those  of  Bates,  et  al.  with  the  use 
of  the  essential  simplifications  that  the  transitions  of  greatest  impor¬ 
tance  are  those  across  the  critical  state  and  that  the  critical  state  can 
be  located  by  minimizing  the  recombination  rate  with  respect  to  the 
critical  quantum  number.  Because  excited- state  population  decreases 
with  increasing  quantum  number,  three -body  recombination  into  an  ex¬ 
cited  state  increases  with  increasing  quantum  number,  collisional  de¬ 
excitation  of  a  given  state  increases  with  increasing  quantum  number, 
and  radiative  recombination  into  a  given  state  decreases  with  increasing 
quantum  number;  they  were  able  to  show  that  the  critical  quantum  number 
is  nearly  equal  to  the  larger  of  the  quantum  numbers  corresponding  to  the 
following  states;  (1)  the  state  at  which  the  collisional  de-excitation  rate 
is  a  minimum  and  (2)  the  state  at  which  collisional  and  radiative  de¬ 
excitation  rates  are  equal. 

Byron,  Bortz,  and  Russel  (Ref.  12)  applied  the  same  philosophy  to 
computation  of  collisional  recombination  rates  for  potassium,  comparing 
the  results  with  those  obtained  for  hydrogen  and  argon  (the  latter  having 
been  computed  in  Ref.  13  from  the  von  Engle  ionization  theory  with  the 
first  excitation  potential  used  instead  of  the  ionization  potential).  Since 
the  critical  state  could  not  be  determined  analytically,  the  following  pro¬ 
cedure  was  used:  at  a  given  electron  temperature  and  an  assumed 
critical  energy  level,  the-energy  levels  and  degeneracies  for  all  levels 
between  the  assumed  critical  level  and  the  level  one  principal  quantum 
number  below  it  were  obtained  from  the  spectral  tables  of  Ref.  14.  The 
transition  rate  across  this  critical  interval  was  then  computed  at  the  given 
electron  temperature.  The  minimum  of  the  curve  of  transition  rates  ob¬ 
tained  by  this  procedure  then  corresponds  to  the  true  critical  energy  level. 
This  tedious  calculation  must  be  repeated  for  each  separate  electron 
temperature. 

The  author  (Ref.  15)  improved  this  procedure  in  the  case  of  alkali 
metals  by  finding  an  analytic  approximation  for  the  degeneracy  averaged 
over  the  substates  within  the  critical  interval,  making  possible  an  analytic 
determination  of  the  critical  energy  level.  The  results  indicated  that  the 
behavior  (with  respect  to  electron  temperature)  of  the  critical  energy  level 
can  be  determined  to  sufficient  accuracy  with  knowledge  required  of  only 
the  ionization  potential,  effective  ground  state  quantum  number,  and  the 
maximum  value  of  the  angulai'  momentum  quantum  number  in  the  interval 
between  the  ground  state  and  the  state  one  principal  quantum  number  above 
ground.  This  procedure  ultimately  yields  ionization  and  recombination 
rates  for  any  element  whose  quantum  mechanical  description  can  be  col¬ 
lapsed  to  that  of  a  one- electron  atom  (whose  potential  deviates  from 
coulombic  form  only  in  the  first  order). 
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Section  II  of  this  report  presents  the  formal  derivation  of  the 
reaction  rates  in  terms  of  parameters  determined  in  Sections  III  and  IV. 

Section  III  is  concerned  with  how  the  average  degeneracy  of  states 
within  the  critical  interval  {the  interval  between  the  critical  energy  level 
and  the  level  one  principal  quantum  number  below  it)  is  related  to  the 
critical  quantum  number.  The  quantum  mechanics  of  a  one- electron 
atom  with  a  potential  having  an  inverse  r*  noncoulombic  term  is  formally 
presented  as  suggested  in  Fermi's  notes  (Ref.  16).  The  form  of  the 
quantum  defect  derived  by  this  procedure  permits  the  degeneracy  to  be 
averaged  over  the  substates  within  the  critical  interval  and  the  result  to 
be  related  to  the  critical  quantum  number  and  the  previously  mentioned 
known  quantities.  In  Section  IV  the  results  of  Section  III  are  used  in  the 
determination  of  the  critical  state:  the  transition  equations  are  mini¬ 
mized  with  respect  to  the  critical  quantum  number,  making  use  of  the 
degeneracy  relation  derived  in  Section  HI.  All  these  parameters  are 
ultimately  used  to  obtain  the  collisional  ionization  and  recombination  rates 
in  Section  V  by  means  of  the  equations  derived  in  Section  II,  and  the  re¬ 
sults  are  added  to  available  radiative  transition  information  in  Section  VI 
to  provide  a  fairly  complete  treatment  of  ionization  and  recombination  in 
moderately  dense  {10l2  _  iol5  electrons/cc)  alkali  plasmas. 

After  completion  of  the  research  reported  herein,  the  author  became 
aware  of  a  calculation  of  cesium  recombination  rates  by  Dugan  (Ref.  17). 

The  Byron  procedure  was  used  together  with  the  first  70  energy  levels  of 
cesium  and  results  obtained  by  use  of  an  electronic  computer  for  electron 
temperatures  ranging  to  10,  OGO^K.  In  the  temperature  range  appropriate 
to  the  physical  restrictions  imposed  by  the  present  author's  approxima¬ 
tion  technique  (0  -  4000°K),  Dugan's  calculations  agree  very  closely  with 
those  of  the  present  theory. 

In  addition  to  the  experiments  of  Refs.  5  and  18,  which  have  dem¬ 
onstrated  the  validity  of  the  general  model  of  collisional  radiative  re¬ 
combination  for  hydrogen  and  helium  afterglows  and  low  density  plasma 
jets,  there  exists  a  limited  amount  of  experimental  data  on  cesium  re¬ 
combination  (Refs.  19  and  20).  In  both  these  references  the  results, 
corrected  for  ambipolar  diffusion  and  other  empirical  complications, 
were  favorably  compared  with  the  theory  of  Ref.  10  for  hydrogen.  The 
present  theoretical  results  for  cesium  and  those  of  Dugan  lie  quite  close 
to  the  aforementioned  calculation  (and,  hence,  the  experimental  data); 
whereas,  the  present  results  for  hydrogen  are  much  less.  From  these 
and  other  considerations  it  is  ultimately  concluded  that  the  present  theory  - 
while  not  treating  collisional  -  radiative  interaction  in  hydrogen  adequately  - 
is  reliable  within  50  percent  (well  within  the  resolution  of  current  experi¬ 
ments)  for  conditions  typical  of  moderately  dense  alkali  plasmas. 
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SECTION  II 

THE  GENERAL  MODEL 


In  this  section  some  properties  of  the  general  model  for  collision- 
dominated  ionization  and  recombination  are  developed  from  an  empirical 
point  of  view  and  incorporated  into  the  mathematics  describing  the  pro¬ 
cess.  Except  for  some  divergence  of  point  of  view  in  the  detailed  char¬ 
acteristics  of  various  theories,  the  derivation  presented  is  more  or  less 
common  to  all  considerations  of  this  problem. 

The  afterglow  experiments  of  Refs.  5  and  6  and  the  spectrographic 
studies  of  low  density  arc  jets  of  Ref,  18  have  demonstrated  the  following 
characteristics  for  hydrogen  and  helium  plasmas;  above  some  critical 
state  all  the  remaining  states  obey  Boltzmann's  distribution  law  --  i.  e. , 
excited- state  densities  (from  spectral  intensities)  are  simple  exponentials 
of  principal  quantum  number.  When  these  population  densities  are  plotted 
on  a  semilogarithmic  graph  against  quantum  number,  the  resulting  slopes 
determine  the  temperahTfe'olf  the  excited- state  distribution;  the  result 
ag]reeT~w]Ih  fre'^~electron  temperatures  determined  from  conductivity^eas- 
ur^imeats.  As  recombination  progresses,  these  states  remain  in  transient 
thermal  equilibrium  with  a  time-varying  electron  temperature.  In  Refs.  10 
and  11  evidence  is  presented  that  the  deceptive  time  behavior  of  the  electron 
concentration,  w'hich  had  for  so  long  prevented  identification  of  three-body 
decay  processes,  can  be  accounted  for  by  the  time  variation  of  the  recom¬ 
bination  rate  resulting  from  the  declining  electron  temperature.  The  net 
result  of  this  dual  relaxation  is  an  appearance  as  two-body  recombination. 

Because  the  various  experiments  cited  previously  have  shown  that  the  ^ 
experimental  excited-state  densities  (for  states  above  some  critical  level)  ^  ? 
taken  together  with  the  experimental  electron  temperatures,  when  substi¬ 
tuted  into  the  Saha  equation,  predict  electron  densities  in  agreement  with 
independent  determinations  of  the  same,  all  theories  of  collision-dominated 
ionization  and  recombination  start  with  the  supposition  that  only  states  be¬ 
low  the  critical  level  need  be  considered  in  determining  the  ionization  and 
recombination  rates.  In  principle  a  finite  matrix  of  transition  equations 
(cut  off  at  the  critical  state)  could  be  solved  simultaneously  for  the  time 
behavior  of  the  excited- state  populations;  however,  it  is  expedient  and 
justifiable  to  consider  transitions  across  all  states  except  one  to  be  very 
fast  in  comparison  to  excitation  or  de- excitation  of  that  one  state.  On  this 
basis  Bates,  et  al.  (Refs.  7  and  8)  computed  the  rate  of  populating  the 
ground  state  by  collisional  and  radiative  transitions.  Byron,  et  al.  (Ref.  11) 
have  taken  the  point  of  view  that  the  critical  transition  is  across  the  interval 
between  the  critical  state  and  the  state  one  principal  quantum  number  below 
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it  (referred  to  herein  as  the  critical  interval).  After  radiative  effects 
are  propely  accounted  for,  their  results  at  high  temperature  are 
practically  identical  with  those  of  Ref.  7.  The  calculation  of  Ref.  10 
is  somewhat  different  on  the  surface  --  classical  ionization  rates  are 
summed  for  all  levels  below  the  critical  level;  however,  since  the  col- 
lisional  transition  cross  sections  automatically  account  for  all  transi¬ 
tions  across  intervals  greater  than  or  equal  to  the  given  interval,  the 
various  models  are  in  reality  equivalent.  Differences  in  the  results  of 
these  investigators  occur,  primarily,  because  of  the  assumptions  made 
in  determining  the  critical  energy  level  (particularly  the  manner  in  which 
it  is  influenced  by  radiative  transitions).  Lesser  differences  are  related 
to  differences  in  cross  sections  used  and  approximations  made  in  eval¬ 
uating  necessary  integrals  involving  the  cross  sections. 


In  formulating  the  equations  of  the  general  model,  the  present  theory 
follows  closely  the  formulation  of  Refs.  11  and  12.  The  point  of  de¬ 
parture  from  the  Byron  model  and  the  primary  contribution  of  the  present 
theory  is  the  development  of  an  explicit  analytical  procedure  for  determin¬ 
ing  the  critical  energy  level. 


The  cross  section  for  exciting  the  "jth"  state  from  some  other  "ith” 
state  by  electron  collisions  is  obtained  by  averaging  the  individual  cross 
sections  over  an  assumed  Maxwell  -  Boltzmann  distribution: 


Q  (f,  “  «i)  = 


fQ  (q  -gt)  Ve 

/ve  f{?e)  d’v. 


(1) 


where  f(v5)  is  the  Maxwellian  velocity  distribution  function,  is  the 
electron  speed  and  Q(ej  -  e;)  represents  the  probability  that  an  electron 
of  speed  will  cause  a  transition  from  energy  level  to  ej.  In  this 
equation  the  thermal  speed  of  the  electrons  has  been  factored  out  of  the 
averaged  cross  section. 


The  state  notation  used  is  in  the  fashion  of  one- electron  atoms 
referenced  to  ground  state  so  that  =  0  refers  to  the  ground  state 
energy  (not  to  be  confused  with  the  permittivity  of  free  space),  fi  to  the 
energy  level  one  principal  quantum  number  above  ground,  etc.  Because 
the  individual  excitation  cross  sections  which  will  be  used  later  are 
semiclassical,  it  must  be  assumed  that  transitions  from  all  substates  of 
a  given  level  to  some  common  upper  state  are  equally  probable,  except 
for  the  weighting  effect  of  the  degeneracy  of  each  substate:  for  example, 
the  collisional  transition  4P  -  5S  is  twice  as  probable  as  the  transition 
4S  -  5S  because  the  P  states  are  twice  as  degenerate  as  the  S  states. 
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In  terms  of  the  averaged  cross  section,  the  rate  of  exciting  the  "jth” 
state  from  all  lower  states  is  written  as 


dn. 


dt 


_  >=)-!  _ 

He  Ve  ^  Hi  Q  (f]  “  fi) 
1=0 


(2) 


where  "i  denotes  number  density  of  the  "ith"  state,  denotes  electron 
number  density,  and  denotes  electron  thermal  speed. 


For  the  computations  to  come  later,  it  will  be  necessary  to  derive 
the  rate  of  crossing  a  given  energy  level.  The  convention  is  now 
adopted  that  all  cross  sections  used  henceforth  represent  the  probability 
of  exciting  transitions  across  an  interval  greater  than  or  equal  to  the 
given  transition  interval.  The  rate  of  such  transitions  is  related  to  the 
cross  sections  as  in  Eq.  (2),  it  being  understood  that  the  notation  refers 
to  crossing  the  "jth''  energy  level,  not  just  to  exciting  the  "jth''  state. 


dnj 

dt 


_  i=)-l 


^  ni  0  (fi  -  fi) 

1=0 


(3) 


dn,  -• 

— —  means  the  rate  of  exciting  the  "jth"  state  and  all  higher  states. 


In  order  to  evaluate  n;  we  must  assume  a  Boltzmann  distribution 
below  a  certain  critical  state.  This  assumption  is  not  entirely  in  accord 
with  experiment.  References  10  and  18  show  that  above  the  critical  state 
Boltzmann  equilibrium  between  actual  states  does  prevail,  the  tempera¬ 
ture  of  the  distribution  being  that  of  the  free  electrons  as  determined  by  -  ^ 
use  of  the  Saha  equation  with  the  measured  free  electron  and  excited- 
state  number  densities.  Below  the  critical  state  the  various  substates  of 
a  given  state  exist  in  equilibrium  at  a  temperature  between  the  gas 
temperature  and  the  free  electron  temperature.  The  states  themselves 
are  not  shown  to  be  in  equilibrium  by  the  spectrographic  data.  Insofar  as 
the  effect  of  this  complication  of  the  basic  model  is  concerned,  we  shall 
accept,  here,  the  alternatives  stated  in  Ref.  7:  either  (1)  a  uniform  state 
distribution  is  maintained  among  states  below  the  critical  energy  level  at 
high  plasma  density  or  (2)  for  a  low  density  plasma  in  which  radiation 
dominates  the  problem  the  state  distribution  is  immaterial.  The  basic 
experimental  character  of  the  distribution  (enhancement  of  intermediate 
state  populations)  is  preserved  in  the  model.  Figure  Vindicates  that, 
depending  on  the  electron  temperature  and  density,  a  state  in  Boltzmann 
equilibrium  with  the  ground  state  is,  generally,  more  densely  populated 
than  when  the  same  state  (as  the  critical  state)  is  evaluated  in  Saha 
equilibrium  with  the  free  electrons. 


A  schematic  diagram  of  the  present  model  is  given  in  Fig.  2.  In  the 
equations  below,  subscript  c  denotes  the  critical  state  and  c  ~  1  refers 
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to  the  state  one  principal  quantum  number  below  the  critical  state.  From 
the  critical  energy  level  on  up  to  the  continuum,  all  states  are  in  Saha 
equilibrium  with  the  free  electrons;  whereas,  all  states  from  ground  state 
to  the  c  -  1  state  are  in  Boltzmann  equilibrium  at  the  free  electron  tem¬ 
perature.  The  energy  difference  between  the  states  c  and  c  -  1  defines  the 
critical  interval  Ae  .  In  terms  of  these  quantities,  the  rate  of  upward 
crossing  the  critical  level  is  written  below,  after  substitution  of  the  Boltz¬ 
mann  distribution  for  states  below  the  critical  state. 


^  1=0  —  1 

He  Vc  ^ 

1=0 


0  ( tc 


<i) 


gi 

Bc-i 


e 


Af 

kT,  ^  kTe 


(4) 


where  g;  is  degeneracy  of  the  "ith"  state,  k  is  Boltzmann's  constant,  and 
Tg  is  the  electron  temperature. 


If  there  are  no  opposing  transitions,  the  rate  of  ionization  is  just  this 
rate  of  upward  transitions  traversing  the  critical  interval.  (The  neglect 
of  opposing  transitions  will  be  examined  further  in  Section  V.  )  Byron's 
assumption  now  will  be  made;  the  difference  between  adjacent  energy 
levels  is  sufficiently  small  in  comparison  to  the  energy  level  (measured 
with  respect  to  ground)  that  the  sum  in  Eq.  (2)  can  be  replaced  by  an 
integral,  the  independent  variable  having  been  transformed  from  state 
index  to  energy  difference  between  (assumed  continuous)  states.  The  re¬ 
sult  is  a  general  statement  of  collision-induced  ionization. 


where  subscript  o  denotes  ground  state  and  g  is  the  average  degeneracy 
of  the  substates  within  the  critical  interval  Ac. 


The  collisional  recombination  expression  can  be  obtained  by  use  of 
detailed  balance,  the  equilibrium  relation  being  the  Saha  equation  (evalu¬ 
ated  at  the  free  electron  temperature). 


(6) 

where  is  the  electron  mass,  h  is  Planck's  constant,  and  <  is  the 
ionization  potential. 


Both  Eqs.  (5)  and  (6)  can  be  written  as  a  common  factor  ?5  multiplied 
by  the  ionization  or  recombination  rate  given  by  the  classical  Thomson- 
Bohr  theory  (see  Ref.  4).  The  classical  theory  considers  a  bound  state 
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electron  is  impulsively  excited  all  the  way  to  ionization  by  collision  with 
a  free  electron  of  appropriate  energy.  Its  results  can  be  expressed  as 
shown  below: 


(—) 

\  ^  *  yio  n  12  a  tio  n 
classiGul 


^  € 

“■  k  T 

ng  no  v'c  e  *  <7^^ 


(7) 


where 


is  the  electronic  charge,  and  is  the  permittivity  of  free  space. 


The  ratio  of  the  sequential  collision- induced  processes  to  the  classical 
impulsive  processes  is 


/in,\ 

^  <lt  j  class  1C  al 


So 


e 

^cl  At 


4Qre,  JE 


(9) 


As  will  be  apparent  later,  S  varies  through  many  orders  of  magnitude  as 
the  electron  temperature  ranges  between  500  and  3000®K.  The  computa¬ 
tion  of  ionization  and  recombination  rates  from  these  formal  equations  de¬ 
pends  on  the  determination  of  the  critical  energy  level  and  the  average 
degeneracy  of  the  substates  within  the  critical  interval.  These  problems 
are  considered  in  detail  in  the  next  two  sections.  The  method  developed 
there  is,  strictly  speaking,  valid  only  in  the  collision-dominated  limit, 
as  it  ignores  the  effect  of  radiative  interactions  in  determining  the  critical 
energy  level.  The  equations  previously  derived,  however,  are  quite 
generally  valid  when  proper  information  is  supplied  concerning  the  critical 
energy  level. 


SECTION  III 

DEGENERACY  OF  STATES  WITHIN  THE  CRITICAL  INTERVAL 

3-1  QUANTUM  MECHANICS  OF  THE  ALKALI  ATOMS 

Inihis  section  the  quantum  mechanical  behavior  of  alkali  atoms  will 
be  determined  in  the  central  force  field  approximation  using  a  procedure 
outlined  in  Fermi's  published  notes  (Ref.  16). 
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On  account  of  the  penetration  of  the  valence  electron  through  the 
closed  shell  of  inner  core  electrons,  the  alkali  atom  energy  levels  can 
be  approximately  calculated  by  use  of  the  potential  shown  below: 


V(0  =  - 


(10) 


where  r  is  the  radial  coordinate  for  the  alkali  valence  electron.  This 
potential  is  clearly  only  approximate.  Both  the  Thomas- Fermi  statis¬ 
tical  theory  (see  Ref.  21)  and  the  old  quantum  theory  (see  Ref.  22) 
indicate  that  an  inverse  r"  term  occurs  because  of  the  polarization  of 
the  core  electrons  by  the  valence  electron  at  large  values  of  r.  The 
polarization  term  is  replaced,  here,  by  the  inverse  r’  term  because 
the  latter  can  be  absorbed  into  the  centrifugal  portion  of  the  separated 
Schroedinger  equation  for  radial  motion  shown  below.  (See.  for  ex¬ 
ample,  Ref.  2  3.  ) 


a, 


1)  R 


0 


(11) 


where 


y  =  2 


Its. 


% 


is  the  reduced  electron  mass,  and  is  the  first  Bohr  radius  for 
hydrogen  (h  =  In  this  equation,  is  an  effective  angular  momentum 

quantum  number;  it  is  not  an  integer.  The  relation  between  I'  and  i,  the 
ordinary  angular  momentum  quantum  number,  is 

£'{£'+!)=  £(£  +  1)  -  ^  (12)  ■ 

The  term  )3  is  a  constant  related  to  the  charge  and  average  radius  of 
the  inner  electronic  core;  it  need  not  concern  us  further.  The  term  y 
is  the  transformed  wave  number  corresponding  to  the  energy  eigenvalue. 
In  order  to  understand  the  characteristics  of  the  alkali  energy  levels,  let 
us  digress  to  the  hydrogen  atom  calculation  (see  Ref.  2  3). 


In  the  special  case  of  the  hydrogen  atom,  the  Schroedinger  equation 
can  be  transformed  by  the  substitution  x  =  yr  with  new  indices  defined  as 


y 


+  £  and  K  ~  2£  +  1 


The  radial  equation  becomes 


X  -JLEl.  +  2  ^  + 


J  “ 


K-  1 


k"  -  1 
4* 


R  =  0 


(13) 


(14) 
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X  K  —  1 

A  standard  transformation,  R  =  e  ^  x  ^  f  (x)  ,  generates  the  form  of 

the  associated  Laguerre  equation,  whose  solution  is  obtained  by  recog¬ 
nizing  that  the  equation  is  identical  to  the  result  of  differentiating  the 
ordinary  Laguerre  equation  «  times  with  respect  to 


.  d*  f  (x)  ,  d  f  (x) 

dx"  dx 


(k  +  1  ~x)  +  (j  -  k)  f(x)  ==  0. 


(15) 


In  terms  of  the  Laguerre  functions,  the  solution  of  the  radial  equation 
for  hydrogen  is 

X- 1  -k; 

„  2  2  d  Lytx) 

R  =  e  X  - h -  (16) 


dx 


Unfortunately,  this  straightforw'^ard  procedure  cannot  be  simply  genera¬ 
lized  for  the  case  of  the  alkali  potential.  Since  the  differentiation  index 
K  must  be  an  integer,  the  above  procedure  cannot  be  followed  when 
K  =  +  ]  is  not  an  integer.  In  this  case  the  radial  solutions  could  be 

obtained  by  numerical  integration  of  Eq.  (8);  our  interest,  however, 
lies  only  in  the  dependence  of  energy  levels  and  degeneracies  on  i’. 
These  relations  can  be  obtained  without  having  explicit  radial  solutions. 


In  analogy  to  the  treatment  of  the  hydrogen  atom,  let  us  require 
that  j  be  a  positive  integer.  If  the  energy  eigenvalues  were  related  to 
an  integral  quantum  number,  the  latter  would  be  determined  as  the 
principal  quantum  number  N  =  j  -  J?.  Following  Fermi,  we  define  an 
effective  principal  quantum  number 


N'  =  N  +  i  - 


(17) 


whose  deviation  from  the  integral  principal  quantum  number  is  the 
quantum  defect  5^.  In  terms  of  these  quantities,  which  characterize 
the  alkali  atoms,  the  energy  eigenvalues  are  written  below: 


V 


2  aj  N' 


(18) 


where 


and 


N'  -  N  -  di 

Si  -  = 


i) 


Since  the  treatment  of  the  alkali  atoms  used  here  assumes  a  central 
force  for  the  field  in  which  the  valence  electron  moves,  the  simple 
degeneracy  of  a  given  substate  is  2(i+  1) ;  when  the  angular  momentum 
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caused  by  electron  spin  is  added,  the  degeneracy  becomes  2(j2+  1).  Now 
that  the  degeneracies  and  energy  levels  of  the  alkali  atoms  have  been 
related  to  appropriate  quantum  numbers,  defined  in  such  a  way  as  to 
collapse  the  alkali  quantum  mechanics  to  hydrogen-like  description,  the 
average  degeneracy  of  the  substates  in  the  critical  interval  can  be  cal¬ 
culated.  The  derivation  is  given  in  the  next  section. 


3.2  AVERAGE  DEGENERACY  APPROXIMATION 


In  this  section  the  average  degeneracy  for  states  within  the  critical 
interval  is  computed  to  the  first  order  in  terms  of  •  By  means  of  the 
approximation  scheme  of  this  section,  the  degeneracy  and  subsequently, 
the  critical  energy  level,  are  related  to  known  parameters  characterizing 
the  alkali  atoms.  The  order  of  approximation  introduced  here  will  be 
followed  in  determining  the  final  ionization  and  recombination  rates,  so 
that  the  present  theory  is  formally  restricted  to  first  order  terms  in  the 
reciprocal  of  the  critical  quantum  number. 


The  derivative  of  the  quantum  defect  is  obtained  from  Eq.  (18): 

d  Sf 

i‘l  +  I  +  — 

2 


The  simple  degeneracies  are  now  averaged  over  the  energy  levels  corre¬ 
sponding  to  substates  within  the  critical  interval  to  obtain  the  average 
degeneracy  I.  In  terms  of  the  integral  below,  g  can  be  written  as 


(19) 


In  this  derivation  denotes  the  maximum  value  of  £.  in  the  critical 
interval.  It  should  be  noted  that  the  Integration  begins  at  some  initial 
value  of  i  with  the  principal  quantum  number  equal  to  -  1 .  As  i  in¬ 
creases  to  the  effective  principal  quantum  number  is  considered  to 
change  from  N'  -  1  to  N'  ;  whether  an  S  state  is  crossed  is  immaterial 
because  N  -  ~  N  -  So  -t-  1 .  The  last  term  in  Eq.  (19)  is  evaluated  by 

use  of  Eq.  (18),  yielding 


/ 


j{  =  0 


Si  (N'  -  1)' 


(20) 
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The  integral  term  can  be  shown  by  use  of  a  second  order  Maclaurin  ex¬ 
pansion  of  S,  to  be  a  second  order  term  in  -W-  With  only  first  order 

N  c 

terms  retained,  the  degeneracy  assumes  the  form  shown  in  Eq.  (21). 


where 


-  1) 

(2  N/-  1  ) 


(21) 


For  comparison,  the  corresponding  degeneracy 
as  Eq.  (22): 

gM  =  (Nc  -  if 


(22) 


-  2  for  hydrogen-like  states. 


1  +  N  ' 


a  varies  from  one  in  the 
for  the  lowermost  states. 


Since  ^ N 

limit  of  very  high  quantum  number  to 

1  +  i  y 

It  will  be  apparent  in  the  next  chapter,  however,  that  the  critical  energy 
level  is  quite  insensitive  to  the  form  of  the  degeneracy  in  the  high  quantum 
number  (low  temperature)  limit.  Because  of  this  fact,  the  alkali  atom  be¬ 
havior  can  be  adequately  characterized,  in  the  present  theory,  by  the  low 
quantum  number  limit;  henceforth,  a  will  be  taken  to  be  its  lower  limit 

1  +  N' 

value,  a  =  ■-  -  - — ,  where  N'  refers  to  ground  state,  effective,  principal 
quantum  number. 


SECTION  IV 

LOCATION  OF  THE  CRITICAL  STATE 


*As  Byron,  et  al.  (Ref.  11)  have  indicated,  the  critical  state  is 
located  so  that  the  rate  of  transition  across  the  critical  state  corresponds 
to  the  slowest  in  a  sequence  of  transitions  from  ground  state  to  the  con¬ 
tinuum.  This  minimal  transition  is  affected  by  both  the  balance  of 
collision-induced  excitation  and  de-excitation  and  the  interaction  of  col- 
lisional  and  radiative  transitions-  In  this  section  the  radiative  effects  are 
ignored  and  the  rate  of  transition  across  the  critical  state  with  respect  to 
the  critical  quantum  number  is  formally  minimized,  making  use  of  the 
parameters  introduced  in  the  preceding  section.  The  general  equation  for 
determining  the  critical  quantum  number  is  shown  below; 


d 


ri  fn  g 


^  in  F  Ik'l’e 

<1  Nf 


1  d  ^  i\  In  hf 

k  '1  e  '1  N  c'  dIV  * 


0 


(23) 
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where 


A( 

<^ol 


oo 


Q(E1 


dE 


(24) 


To  solve  this  equation,  let  us  use  the  Gryzinski  (Ref.  9)  cross  sections 
for  collision- induced  excitation,  as  linearized  by  Byron,  et  al.  (Ref.  12). 
The  various  cross  sections  will  be  considered  in  more  detail  in’ a  later 
section.  The  parameters  to  be  inserted  into  Eq.  (23)  are  listed  below: 


d  in  g  f 

Nc'  \ 

r2iNc'  -  (a+  2)  1 

d  in  N/  \ 

Nc'  -  1  / 

L 

Nc'  -  (c  +  ]  )  J 

No'’ 

d  N,,' 

N/’ 

J 

^  d  in  At  \ 

d  in  F  e/ 

d  In  He  ) 

d  In  Nj  ' 

1 

/  ''Te/ 

d  in  At 

2 

3Nc'(Nc'-  1)  +  1 

d  In  Nc' 

(Ne'  - 

1) 

2Nc'-  1 

(25) 


At  this  stage  it  is  advantageous  to  split  the  critical  quantum  number 
according  to  the  following  transformation: 


kT. 


and  = 


f  -  Ee 


n;  ’  ’ '  -  -  kT. 

The  two  new  parameters  ^  and  C  are  of  direct  physical  interest;  is 
proportional  to  the  shift  of  the  critical  energy  level  down  from  the  ioniza¬ 
tion  potential,  and  is  the  dimensionless  free  electron  mean  energy. 

The  normalizing  factor  in  the  expression  for  is  related  to  the  Rydberg 
constant  for  the  element  in  question.  The  normalization  has  the  effect  of 
removing  any  dependence  on  atomic  mass,  so  that  the  true  effect  of  dif¬ 
ferences  in  atomic  structure  can  be  ascertained.  In  terms  of  these 
quantities  and  to  the  first  order  approximation  with  respect  to  ^  C ,  the 
equation  finally  determining  the  critical  energy  lev.el  is  written  below  as 
a  simple  quadratic; 


=  0 


(26) 


where 


1  d  In 


X  =  — 

2  d  in  N/ 
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d  in  g  6 

d  in  Nc'  ^  2  - 
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The  solution  to  this  equation  is  presented  in  Fig.  3.  The  tempera¬ 
ture  ranges  from  zero  to  about  2500®K.  The  present  theory  apparently 
compares  favorably  (within  about  10  percent  at  high  temperature)  with 
the  specific  determinations  of  Ref.  12.  From  quantities  stated  in  Eq.  (25), 
one  can  see  that  at  high  critical  quantum  numbers,  the  degeneracy  deriva¬ 
tive  is  essentially  hydrogen-like  and  is  independent  of  a.  This  behavior 
accounts  for  the  previous  remark  that  the  lower  limit  of  a  is  sufficient  to 
characterize  the  behavior  of  the  various  one-electron  substances.  As  is 
apparent  in  Fig.  3,  the  theory  and  Byron's  points  converge  to  hydrogen¬ 
like  behavior  at  the  high  quantum  number  (low  temperature)  limit.  It  is 
noteworthy,  however,  that  at  the  upper  end  of  the  curves,  where  the  dis¬ 
crepancies  between  theory  and  Byron's  points  occur,  the  disparity  exists 
for  hydrogen  (whose  degeneracy  is  exact)  as  well  as  for  the  alkalies. 

The  reason  for  this  discrepancy  is  the  order  of  approximation  inherent  in 
Eq,  (26).  In  this  equation  consistency  with  the  order  of  approximation  in 
the  degeneracy  computation  has  been  maintained  because  inclusion  of 
second  order  terms  in  Eq.  (26)  (with  only  first  order  terms  in  the  de¬ 
generacy)  leads  to  a  rapid  divergence  of  C-  Consistent  inclusion  of  second 
order  terms  in  both  the  degeneracy  and  the  determination  of  the  critical 
state  should  alleviate  this  discrepancy;  however,  such  additional  and 
arduous  mathematical  labor  is  not  worthwhile  because  error  in  determina¬ 
tion  of  the  critical  energy  level  propagates  into  the  final  ionization  and 
recombination  rates  only  as  the  square  of  the  relative  error  in  It 
should  be  noted,  however,  that  errors  in  the  degeneracy  propagate  linearly 
into  the  final  rate  expressions.  The  significance  of  these  errors  in  the 
approximations  will  be  more  fully  explored  in  Section  V. 

Before  turning  to  the  final  computation  of  the  ionization  and  recombina¬ 
tion  rates,  it  is  of  some  interest  to  compare  various  estimates  of  the 
critical  energy  level  behavior.  The  present  results  indicate  approximately 

i  N' 

:  whereas,  Hinnov  and  Hirschberg  (Ref.  10)  found = 

This  apparent  discrepancy  can  be  explained  on  the  basis  of  radiative 
transitions  across  the  critical  interval.  Hinnov  and  Hirschberg  determined 
the  critical  state  (for  electron  mean  energy  <  4  ev)  by  equating  forward  and 
reverse  transitions  across  the  critical  interval,  requiring  Boltzmann  equi¬ 
librium  between  the  states  c  and  c-l,  and  using  very  high  critical  quantum 
number  approximations  to  solve  the  resulting  transcendental  equation  for 
N'c  ,  These  assumptions  are  inconsistent  with  those  of  the  present  model, 
wherein  the  critical  state  is  assumed  in  Saha  equilibrium  with  the  free 
electrons  and  the  very  high  quantum  number  approximation  is  not  appro¬ 
priate. 
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The  Hinnov  and  Hirschberg  critical  state  calculation  actually  corre¬ 
sponds  to  assuming  collision-induced  transitions  up  to  and  across  a 
radiation-limited  critical  interval.  It  is  applicable  whenever  radiation 
is  the  predominant  mechanism  of  crossing  the  critical  interval;  this 
predominance  occurs  late  in  the  afterglow  and  corresponds  in  their  ex¬ 
periment  to  electron  mean  energy  <  a ev.  Above  this  energy  the  free 
electron  density  is  sufficient  to  maintain  a  collision-limited  critical  state, 
and  Hinnov  and  Hirschberg  use  a  different  formula  (accounting  for  both 
collision-induced  and  radiative  transitions  between  specified  states)  pro¬ 
posed  in  Ref.  6  and  repeated  in  Ref.  10.  This  calculation  agrees  numeri¬ 
cally  with  the  present  theory  and  their  spectroscopic  data  (2^p-N^D  lines 
in  helium),  yielding  the  critical  quantum  number  =  4.  Below  4®''>  their 
low  temperature  formula  overestimates  and  the  present  theory  under¬ 
estimates  the  critical  quantum  number.  In  this  region  the  critical  quantum 
number  is  best  estimated  by  two  other  theories.  Both  a  simple  calculation 
by  Ref.  11  and  a  much  more  elaborate  one  quoted  in  Ref.  24  give 
y 

lying  between  the  present  results  and  those  of  Ref.  10. 

In  summary  of  these  results,  one  c&n  say  that  the  present  theory  and 
the  numerical  points  of  Ref.  12  underestimate  the  critical  quantum  num¬ 
ber  whenever  radiative  effects  are  important.  Although  such  radiative 
predominance  can  often  occur  in  hydrogen  and  helium  laboratory  plasmas, 
it  is  most  uncommon  in  alkali- seeded  plasmas,  so  that  the  present  theory 
provides  a  good  description  for  plasmas  of  interest  in  practical  MHD 
devices. 

To  date  there  have  been  no  experimental  determinations  of  the  critical 
level  with  the  exception  of  the  previously  mentioned  spectroscopic  studies. 
Observation  of  the  time  variation  of  various  spectral  line  intensities 
coupled  with  measurements  of  electron  temperature  and  density  variation 
and  comparison  with  theoretical  calculations  should  provide  a  more  definite 
determination  of  the  critical  level  behavior,  if  such  is  needed.  This 
possibility  is  considered  in  Section  VII. 

SECTION  V 

FINAL  COLLISIONAL  IONIZATION  AND  RECOMBINATION  RATES 
5.1  GENERAL  THEORY 

Now  that  the  critical  energy  level  behavior  has  been  determined,  the 
collisional  ionization  and  recombination  rates  can  be  determined  from 
Eqs.  (5)  and  (6),  Since  the  factor  S  is  common  to  both  the  ionization  and 
recombination  rates,  only  5  will  be  formally  considered  here.  In  order 
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to  simplify  the  computation  of  5  in  a  fashion  consistent  with  the  deter¬ 
mination  of  the  average  degeneracy  "g,  one  must  determine  the  minimum 
quantum  number  consistent  with  the  first  order  approximation  (in 

terms  of  inherent  in  g.  Higher  order  terms  can  be  eliminated  from 
consideration  in  S  by  using  the  criterion  developed. 

The  basic  inequality  expressing  the  limit  of  the  first  order  degeneracy 
approximation  is 


1 


< 


(27} 


This  inequality  imposes  the  following  restriction  on  the  value  of  the 
critical  quantum  number: 


or 


< 


r  .  .  1 

% 

L  (i  + 

-  f  1  + 

V  4  y 

V  2/ 

L  _ 

Nc'  > 

*  -i-)" 

/ 

(28) 


These  limits  of  the  mathematical  approximations  also  agree  approxi¬ 
mately  with  the  actual  quantum  numbers  for  the  states  one  principal 
quantum  number  above  ground:  for  cesium  the  mathematical  limit  is 
2.  63  and  the  actual  quantum  number  is  2.  93;  for  potassium  the  mathe¬ 
matical  limit  is  3.  12  and  the  actual  quantum  number  is  2.93.  It  is  thus 
apparent  that  for  potassium  the  limit  of  validity  of  the  mathematical 
approximations  is  reached  shortly  before  the  critical  level  has  lowered 
to  the  level  one  principal  quantum  number  above  the  ground  state;  where¬ 
as.  for  cesium,  the  mathematical  approximations  are  adequate  at  the 
physical  limit. 

For  both  these  elements  the  assumptions  of  the  theory  actually  be¬ 
come  invalid  when  the  critical  state  approaches  the  limit  rnentioned. 
Specifically,  the  assumption  of  equal  transition  probability  (except  for 
degeneracy)  for  transitions  originating  from  the  various  substates  and 
crossing  the  critical  state  becomes  quite  questionable  on  account  of  the 
large  differences  in  spacing  of  the  substates  which  constitute  the  interval 
between  the  ground  state  and  the  state  one  principal  quantum  number  above 
ground.  The  procedure  adopted  here  is  to  calculate  the  specific  limit  of 
S  associated  with  the  critical  state  fixed  at  the  state  one  principal  quantum 
number  above  ground.  Two  extremes  of  the  degeneracy  are  assumed 
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corresponding  to  the  assumption  that  (1)  the  only  transition  of  importance 
is  that  originating  at  the  ground  state  and  (2)  transitions  originating  from 
all  states  within  the  interval  are  equally  important.  The  letter  assump¬ 
tion  produces  transition  rates  somewhat  larger  than  the  former.  Pre¬ 
sumably,  if  the  correct  distribution  of  transitions  were  known  (from  a 
quantum  mechanical  calculation  of  transition  probabilities)  the  associated 
transition  rates  should  at  some  temperature  be  tangent  with  an  extrap¬ 
olation  of  the  original  theory  into  regions  where  it  is  in  doubt.  * 


With  all  these  restrictions  in  mind,  let  us  now  write  the  form  of  S 
consistent  with  the  physical  limitations  of  the  model  as  below: 


.-1)1 

f  j  1_  a  \  a/J 
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F 


(29) 


and 


^Limit  ( ) 


UtJ 


(30) 


where  (j  is  the  energy  level  of  the  state  one  principal  quantum  number 
above  ground. 


5.2  COLLISIONAL  EXCITATION  CROSS  SECTIONS 

Presently,  there  exists  no  completely  satisfactory  method  for 
determining  collision- induced  transition  probabilities.  According  to 
Ref.  25,  quantum  theory  is  capable  of  treating  this  problem  exactly; 
however,  the  difficulty  of  the  mathematics  involved  has  prevented  the 
successful  exploitation  of  quantum  mechanical  techniques  for  all  but  the 
simplest  systems.  Even  for  hydrogen  the  Born  approximation  predicts 
low  energy  results  which  overestimate’  experimental  results  by  about 
an  order  of  magnitude;  at  high  energies  the  two  agree  fairly  well,  as 
would  be  expected  according  to  the  nature  of  the  Born  approximation. 
Various  additional  sophistications  (such  as  inclusion  of  exchange  effects 
and  consideration  of  coupling  of  states)  have  lowered  the  cross  section 
somewhat,  at  the  expense  of  such  mathematical  complexity  that  exten¬ 
sion  to  more  complicated  systems  than  hydrogen  is  most  unlikely. 
Omidvar  (Ref.  26)  has  developed  a  procedure  (using  the  electronic  com¬ 
puter)  for  calculating  elastic  and  inelastic  cross  sections  in  hydrogen. 
His  results  lower  the  theoretical  cross  sections  (from  the  Born  value) 
to  values  in  agreement  with  the  experiments  of  Ref.  27  but  considerably 


*The  author  acknowledges  Dr.  L,  Ring  for  the  demonstration  of  this 
point. 
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higher  than  those  of  Ref.  28.  He  indicates  that  inclusion  of  two  more 
states  in  the  coupling  consideration  (he  considered  four)  will  reduce  the 
results  about  30  percent.  Although  these  theoretical  trends  are.  clearly, 
in  the  right  direction,  understanding  of  collision- induced  transitions  is 
far  from  complete  --  one  unexpected  difficulty  being  the  lack  of  agree¬ 
ment  of  the  experiments  mentioned. 

In  an  effort  to  extract  some  information  concerning  quantum  mechani¬ 
cal  peculiarities  of  the  problem  without  actually  solving  it,  Seaton 
(Ref.  -2  5)  has  assumed  the  coiy.sional  excitation  process  involves  emis¬ 
sion  and  reabsorption  of  a  virtual  photon.  From  this  model  he  has  shown 
the  collisional  transition  probability  for  optically  allowed  transitions  to  be 
proportional  to  the  oscillator  strength  for  the  analogous  optical  transition; 
hence,  one  may  expect  some  sort  of  selection  rule  to  govern  collisional 
transitions,  just  as  is  true  for  radiative  transitions. 

In  the  absence  of  such  selection  rules,  one  can  only  use  classical 
cross  sections,  assuming  the  most  probable  transition  is  that  for  which 
AH  =  0  and  partially  accounting  for  the  other  transitions  (across  smaller 
intervals  than  the  primary  one  just  mentioned)  by  assigning  them  the  same 
a  priori  probability,  except  for  the  weighting  effect  of  the  various  de¬ 
generacies.  Still,  one  must  choose  a  cross  section  for  the  lumped  transitions. 

Gryzinski  (Ref.  9)  has  considered  classical  collisions  between  free 
electrons  and  bound  state  electrons  moving  in  allowed  orbits  and  with 
velocities  compatible  with  quantum  mechanical  limitations.  His  cross  sec¬ 
tions  have  been  evaluated  by  Kingston  (Ref,  29)  for  hydrogen  and  com¬ 
pared  with  various  similar  approximations,  experimental  data  (directly 
obtained  from  ionization  experiments  and  indirectly,  by  comparison  of 
collisional- radiative  recombination  theory  with  experiments  in  hydrogen 
and  helium)  and  the  Born  approximation.  He  concludes  that  the  cross  sec¬ 
tions  for  ionization  are  accurate  within  about  20  percent;  whereas,  those 
for  excitation  are  good  within  at  least  a  factor  of  two.  At  high  energy,  the 
classical  theory  underestimates  the  Born  results,  because  the  quantum 
theory  effectively  considers  impact  parameters  extending  to  infinity.  The 
Born  cross  section  falls  off  as  (/nE)  at  large  energies;  the  classical 
cross  section  goes  as  1/E  ,  At  low  energies,  classical  theory  overestimates 
reality  somewhat,,  because  classical  theory  assumes  the  excitation  probabil¬ 
ity  is  one  whenever  the  impact  parameter  of  the  incoming  electron  lies 
within  a  specified  range  (see  Ref.  25).  In  this  respect  the  cross  sections  in 
Ref.  9  are  superior  to  those  of  the  Thomson-Bohr  (see  Ref.  4)  theory  be¬ 
cause  the  averaging  over  all  locations  of  the  bound  electron  has  the  effect  of 
"smearing  out"  the  impact  parameter  in  a  fashion  reminiscent  of  quantum 
mechanics.  (The  Thomson-Bohr  theory  represents  the  limiting  case  of 
Gryzinski's  theory  when  the  energy  of  the  bound  electron  is  zero.  ) 
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Gryzinski  {Ref.  9)  writes  the  cross  section  for  excitation  of  a  bound 
state  electron  with  energy  Ej  through  an  interval  E  by  collision  with  an 

electron  of  energy  —  as 


Q(E)  = 


(31) 


where  gj  is  a  function  of  the  colliding  electron  energy,  the  bound 
electron's  initial  energy,  and  the  transition  interval.  This  function  is 
plotted  in  Fig.  4.  Because  of  the  difficulty  of  integrating  this  function 
over  a  Maxwell-Boltzmann  distribution  of  electron  velocities,  Byron, 
et  al.  (Ref.  12)  use  a  linear  approximation  to  the  first  part  of  the 
Gryzinski  curves,  the  slope  being  a  weak  function  of  the  electron  energy 
on  account  of  the  dependence  of  the  bound  electron's  initial  energy  on  free 
electron  energy.  A  different  procedure  is  used  in  the  present  calcula- 
tionsj  because  the  Thomson-Bohr  cross  sections  can  be  easily  inte¬ 
grated  over  a  Maxwell-Boltzmann  distribution,  these  cross  sections 
are  adjusted  by  a  factor  of  0,  8  to  bring  them  into  approximate  agreement 
with  the  low  energy  portion  of  the  curve  marked  Ej/E  =  i.  The  high 
energy  portion  of  these  curves  is  now  underestimated;  however,  the  ex¬ 
ponential  tail  of  the  Maxwellian  distribution  will  suppress  these  high 
energy  transitions,  so  that  the  effect  should  not  be  noticeable.  This 
treatment  probably  underestimates  the  transition  probabilities  for  high 
quantum  number,  but  certainly  by  less  than  a  factor  of  two.  Indeed  the 
final  results  at  low  temperature  agree  well  with  both  experiments  and 
other  theories,  as  will  be  apparent  shortly. 


In  the  equations  below  the  Thomson-Bohr  cross  sections  and  the 
linearized  Gryzinski  cross  sections  are  stated  and  integrated  over  a 
Maxwell-Boltzmann  distribution; 


where 
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where  B  is  Byron's  dimensionless  slope  of  Gryzinski's  function  gj.-  The 


cross-section  integrals 
cross  sections; 


are  nov/  determined  for  the  preceding 


and 


(34) 


(35) 


All  pertinent  quantities  having  been  determined,  we  can  now  calculate 
the  ionization  and  recombination  rates  from  ^  (determined  by  Eqs.  (26) 
and  (29)  or  (30)  and  the  relations  for  ionization  and  recombination  coef¬ 
ficients  stated  below: 
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In  the  following  section,  recombination  coefficients  computed  by  use 
of  the  modified  Thomson-Bohr  cross  section  are  discussed.  Though  they 
are  not  presented,  results  obtained  from  the  linearized  Gryzinski  cross 
sections  are  almost  identical. 

5.3  COLLISIONAL  RECOMBINATION  RESULTS 

The  final  results  of  the  collision- induced  recombination  study  are 
presented  in  Fig.  5  as  the  variation  with  electron  temperature  of  the 
effective  three-body  recombination  coefficients.  Only  hydrogen,  cesium, 
and  potassium  have  been  considered  because  the  behavior  of  lithium  and 
sodium  is  expected  to  resemble  that  of  potassium;  similarly,  cesium  and 
rubidium  should  recombine  at  about  the  same  rate. 

One  notices  immediately  that  the  various  theories  for  different  ele¬ 
ments  agree  within  the  resolution  of  most  present-day  experiments  at  low 
temperatures:  however,  at  high  temperatures  there  is  considerable 
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divergence.  The  basic  differences  between  the  theories  of  Refs.  11,  12, 
and  the  present  author  are  related  to  the  first  order  approximation  for 
finite  critical  quantum  number  effects  inherent  in  the  present  theory. 

In  turn,  all  these  theories  differ  from  that  of  Ref.  10  by  the  basic  differ¬ 
ences  in  the  computation  of  the  critical  quantum  number.  The  present 
theory,  in  particular,  computes  hydrogen  recombination  rates  at  high 
temperature  an  order  of  magnitude  less  than  the  work  of  Ref.  10. 

The  present  results  agree  well  with  the  theoretical  studies  of  cesium 
reported  by  Dugan  {Ref.  17)  which  came  to  the  author's  attention  after 
completion  of  the  present  study.  Dugan's  calculation  extends  to  10,  000®K. 
Because  his  method  of  determining  the  critical  state  (by  means  of  the 
spectral  data  and  the  use  of  an  electronic  computer)  does  not  exclude  the 
collapse  of  the  critical  interval  into  transitions  which  violate  hJU  =  o, 

Dugan's  calculations  include  contributions  from  critical  states  lower  than 
the  state  one  principal  quantum  number  above  ground.  Consequently,  his 
high  temperature  results  are  somewhat  higher  than  those  of  the  present 
theory.  At  the  other  end  of  the  temperature  scale  Dugan's  results  are, 
generally,  slightly  lower  than  those  of  the  present  theory  because  he  in¬ 
cludes  forward  as  well  as  reverse  transitions  across  the  critical  state. 

In  the  present  theory  the  tendency  to  overestimate  the  net  transition  rate 
because  of  neglect  of  opposing  transitions  is  partially  compensated  by  the 
underestimation  of  the  collision  cross  section  for  transitions  between  high 
quantum  number  states  mentioned  in  the  previous  section.  At  any  rate, 
for  most  of  the  temperature  range  shown  Dugan's  numerical  solutions  and 
the  present  calculations  are  practically  indistinguishable. 

Some  confirmation  of  the  present  results  can  be  obtained  from  the  re¬ 
combination  experiments  of  Refs.  19  and  20  in  near  thermal  cesium 
plasmas  and  Ref.  30  in  potassium- seeded  argon.  The  data  of  Harris 
(Ref.  20)  are  not  shown  in  the  figure  because  Harris  included  an  empirical 
correction  for  ambipolar  diffusion  and  molecular  ion  dissociative  recom¬ 
bination.  He  added  corrections  for  these  effects  to  recombination  mean 
free  times  obtained  from  the  theory  of  Ref.  10;  the  results  agreed  well  with 
the  experimental  measurements  for  electron  densities  greater  than 
5  X  10^^/cc.  Below  this  density  he  found  faster  recombination  than  that 
predicted  by  theory.  In  these  experiments  the  electron  temperature  range 
was  from  1250  to  1700‘’C.  In  this  range  the  present  theory  for  cesium 
agrees  with  the  calculation  of  Ref.  10  for  hydrogen  within  about 
20  percent--far  within  the  experimental  inaccuracy.  This  agreement  with 
the  experimental  data  of  Ref.  20  is  regarded  as  a  sensitive  test  of  the 
present  method  of  determining  the  critical  energy  level  at  electron  tem¬ 
peratures  less  than  2000°K. 
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At  high  temperatures,  v/here  the  finite  quantum  number  effects  be¬ 
come  appreciable,  the  present  theory  for  cesium  agrees  fairly  well  with 
the  data  of  Ref,  ID;  however,  the  density  dependence  (resulting  from 
radiative  effects)  of  these  experiments  produces  sufficient  scatter  in  the 
data  that  conclusive  comparison  is  not  possible  at  this  stage.  Radiative 
effects  in  a  collision-dominated  plasma  will  be  accounted  for  in  the  next 
section,  and  further  comparison  of  experiment  and  theory  must  await 
development  of  these  effects. 

The  data  of  Ref.  30  (which  is  quite  a  recent  publication)  show  con¬ 
siderable  scatter;  however,  the  results  generally  confirm  the  theoreti¬ 
cal  predictions.  The  numerical  computations  of  Byron  et  al.  (Ref.  12) 
more  nearly  fit  the  average  of  the  data  than  do  the  present  results  be¬ 
cause  of  inaccuracies  in  the  degeneracy  approximation  to  be  discussed 
below.  The  theoretical  uncertainties,  nevertheless,  are  well  within  the 
experimental  scatter. 

As  the  critical  state  approaches  the  state  one  principal  quantum 
number  above  ground,  the  recombination  coefficient  assumes  the 
limiting  form  discussed  earlier.  Such  a  transition  is  quite  evident  in 
the  cesium  curve.  Extrapolation  of  the  original  curve  provides  a 
reasonable  estimate  between  the  two  lower  limit  curves  shown.  In  the 
case  of  potassium,  the  mathematical  assumptions  in  the  degeneracy 
approximation  fail  before  the  physical  limit  is  reached,  and  the  high 
temperature  portion  of  the  curve  underestimates  the  lower  limit  con¬ 
siderably  on  account  of  the  progressive  deterioration  of  the  degeneracy 
approximation.  It  can  be  shown  that  the  smaller  a  is,  the  better  is  the 
approximation  for  g-  Since  cesium  includes  the  5D  states  in  the  interval 
between  ground  state  and  the  state  one  principal  quantum  number  above 
ground  and  the  potassium  lowermost  critical  interval  has  no  such  II 
states,  Ocesium  <  “potassium’  hence  the  degeneracy  approximation  for  cesium 
remains  good  after  the  potassium  approximation  has  failed.  Physically, 
the  presence  of  the  D  states  also  means  that  the  lower  limit  curve  for 
cesium  marked  'all  transitions'"  is  more  appropriate  than  the  lower  limit 
curve  marked  "primary  transition’’  in  describing  the  high  temperature 
behavior  of  cesium  recombination.  The  trend  of  the  extrapolated  curve 
agrees  with  this  observation. 


SECTION  VI 

RADIATIVE  CONSIDERATIONS  AND  FINAL  COMPARISON 
BETWEEN  EXPERIMENTS  AND  THEORY 


Although  the  problems  of  interest  in  this  study  are  primarily  con¬ 
cerned  with  collisional  effects,  it  is  necessary  to  consider  the  implications 
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of  radiative -collisional  interaction  over  a  range  of  optical  opacities 
before  the  final  comparison  of  theory  and  experiment  can  be  made.  In 
addition  to  the  small  direct  effect  of  radiative  recombination  (or  photo¬ 
ionization),  the  ionization  and  recombination  rates  are  affected  by  the 
influence  of  radiation  on  the  distribution  of  excited  states  (hence,  the 
location  of  the  critical  state).  In  addition,  one  may  question  the  effect 
of  additional  cascading  between  lower  states  resulting  from  radiation. 

(This  effect  has  been  taken  into  account  in  the  second  calculation  for 
hydrogen  at  high  electron  temperatures  in  Ref.  10.) 

The  assumption  of  Boltzmann  equilibrium  between  states  lower 
than  the  critical  state  is  equivalent  to  an  infinite  rate  of  cascading  (up 
and  down)  so  that,  within  the  context  of  the  present  theory,  the  effect 
of  line  radiation  can  be  considered  only  to  the  extent  of  its  effect  on  the 
location  of  the  critical  state  and  the  enhancement  of  transitions  across  - 
the  critical  interval.  It  should  be  noted  that  if  the  plasma  were  optically 
thick  to  all  line  radiation,  a  Boltzmann  distribution  of  lower  states  would 
be  assured  on  account  of  the  mutual  requirements  implied  by  blackbody 
radiationj  the  balance  between  absorption,  induced  emission,  and 
spontaneous  emission;  and  thermal  equilibrium  between  excited  states 
and  the  radiation  field.  Similarly,  in  the  limit  of  collision  domination, 
thermal  equilibrium  between  free  electrons  and  valence  electrons  implies 
a  quasi- equilibrium  distribution  of  states,  if  the  free  electrons  undergo 
sufficient  elastic  collisions  to  maintain  a  Maxwell-Boltzmann  distribution 
of  free  electron  speeds.  (This  point  is  mentioned  in  Refs.  7  and  8  and  in 
Ref.  31  in  a  discussion  of  local  thermodynamic  equilibrium.  ) 

The  effect  of  line  radiation  in  determining  the  distribution  of  states 
and  the  critical  state  location  can  be  assessed  by  comparing  radiative  and 
collisional  transition  probabilities  across  the  critical  interval.  Since  the 
collisional  transitions  considered  in  this  report  correspond  to  optically 
forbidden  single  transitions,  the  competition  between  collisions,  two  stage 
allowed  radiative  transitions  across  the  same  interval  as  the  collisions, 
and  optically  allowed  single  transitions  across  some  interval  nearly  equal 
to  the  collisional  interval  must  be  considered.  Since  the  line  intensities 
are  usually  represented  by  the  Einstein  transition  probabilities  per  unit 
time,  the  relative  effectiveness  of  radiation  versus  collisions  can  be 
represented  as  the  quotient  of  the  radiative  transition  probability  per  unit 
time  and  the  collisional  transition  probability  per  unit  time.  For  the  sake 
of  convenience,  it  will  be  assumed  that  all  transitions  crossing  a  given 
critical  interval  are  equally  energetic  but  differ  in  probability  on  account 
of  having  different  oscillator  strengths.  The  composite  oscillator  strength 
for  a  dual  transition  will  be  defined  as  the  square  root  of  the  product  of 
oscillator  strengths  for  the  individual  transitions.  The  oscillator  strength 
for  all  transitions  across  a  given  critical  interval  will  be  computed  by 
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summing  the  oscillator  strengths  for  all  complete  transitions  across 
that  interval.  The  necessary  mathematical  relations  for  computing 
these  effects  are  written  below. 


The  Einstein  transition  probability  coefficient  for  emissive  transi¬ 
tions  across  the  critical  interval  is  obtained  from  Ref.  32  in  Kramers' 
form  (the  Gaunt  factor  is  taken  to  be  one): 
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where  is  the  frequency  of  the  transitions  across  the  critical 
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interval,  c  is  the  speed  of  light,  and  represents  the  composite 

oscillator  strength  for  all  possibilities  of  emissive  crossing  the  critical 
interval.  The  analagous  collisional  transition  probability  (per  unit  time) 
is 
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The  ratio  of  radiative  to  collisional  transitions  is  obtained  from 
Eqs.  (38)  and  (39)  after  use  of  Boltzmann's  and  Saha's  relations: 
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It  is  the  same  for  both  excitation  and  de- excitation. 

Table  I  shows  the  compilation  of  oscillator  strengths  for  crossing 
various  critical  intervals.  These  results  presented  in  Ref.  33  for  cesium 
are  only  approximate  and  caution  must  be  exercised  in  applying  them  be¬ 
cause  the  Thomas-Kuhn  sum  rule  (see  Ref.  21)  is  occasionally  violated 
unavoidably.  The  oscillator  strengths  whose  values  are  listed  as  1+  were 
quoted  as  greater  than  one  in  the  original  article. 

Table  II  shows  the  values  of  -Y  obtained  from  the  quoted  oscillator 
strengths  and  Eq.  (40),  It  is  apparent  that  radiative  transitions  outweigh‘d 
collisional  crossing  of  the  interval  between  ground  state  and  the  state  one 
principal  quantum  number  above  ground;  w'hereas.  collisional  transitions 
are  more  important  than  radiation  in  crossing  higher  critical  intervals. 
Calculations  in  Ref.  34  indicate  that  for  transitions  from  excited  states  to 
ground  most  laboratory  plasma  are  optically  thick,  except  for  the  (Doppler 
and  pressure  broadened)  wdngs  of  the  resonance  lines.  On  the  basis  of 
those  results  it  will  be  assumed  in  the  present  computation  that  excited 
state  to  ground  transitions  are  completely  absorbed  and  that  other  line 
transitions  represent  small  perturbations  of  the  collisional  transition  rate. 
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The  enhancement  of  recombination  by  line  radiation  is  ultimately 
computed  by  multiplying  the  collisional  results  by  1  +  X-  When  these 
computations  are  compared  with  the  experimental  data  of  Ref.  19  an  un¬ 
expected  difficulty  is  encountered:  at  about  2300°K  the  critical  state  is 
8  S;  the  next  higher  critical  state  8  P  occurs  at  about  1900°K.  Since  the 
range  of  the  data  of  Ref.  19  is  2000“K  5  1  2320“K,  it  is  evident  that 

some  mixing  of  critical  states  must  occur  at  the  intermediate  tempera¬ 
tures.  Equal  sharing  would  mean  that  2100®K  should  be  the  dividing  line 
between  temperatures  associated  with  the  two  critical  states.  Since  the 
transitions  associated  with  the  8  P  critical  state  are  about  three  times 
more  probable  than  those  across  the  8  S  critical  state  at  their  respective 
temperatures,  the  line  of  demarcation  will  arbitrarily  be  placed  at 
2200'’K.  For  temperatures  above  2200‘^K,  the  critical  state  is  taken  to  be 
8  S,  for  temperatures  below  22  00^K,  the  critical  state  is  taken  to  be  8  P  . 
Although  this  oversimplification  merely  specifies  the  range  over  which 
each  of  the  two  critical  states  predominates,  rather  than  considering  the 
mixing  of  critical  states  at  a  given  temperature,  the  procedure  does  pro¬ 
vide  a  guideline  for  estimating  radiative  enhancement  for  cesium  recom¬ 
bination  rates. 

Figure  6  shows  the  final  comparison  of  the  present  theory  with  the 
available  experiments.  The  straight  lines  represent  collisional  theory 
without  radiative  effects.  In  the  radiative  computation  direct  radiative 
recombination  has  been  ignored,  as  it  is  nearly  constant  over  the  tem¬ 
perature  range  and  on  the  order  of  10"  12  cm^/sec.  Considering  the 
scatter  of  the  data,  the  agreement  between  theory  and  experiment  is 
satisfactory;  moreover,  within  the  uncertainties  present,  purely  colli¬ 
sional  theory  provides  a  good  description  of  the  recombination  process. 


SECTION  VII 
CONCLUDING  REMARKS 


This  report  has  presented  a  theoretical  study  of  ionization  and  re¬ 
combination  in  alkali  plasmas  having  electron  densities  of  10l2  -  IQl^ 
electrons/cc.  The  transition  model  used  assumes  the  existence  of  a 
critical  state  in  equilibrium  with  free  electrons.  The  net  ionization  or 
recombination  rate  is  equal  to  the  transition  rate  across  the  critical 
interval  defined  by  the  critical  energy  level  and  the  level  one  principal 
quantum  number  below  it. 

The  primary  contribution  of  this  investigation  has  been  the  determina¬ 
tion  of  the  critical  energy  level  behavior  (with  respect  to  electron  tempera¬ 
ture).  From  an  approximate  quantum  mechanical  description  of  the  alkali 
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atom  as  one  electron  revolving  in  a  (first  order)  noncoulombic  central 
force  field,  the  energy  level,  quantum  defect,  and  simple  degeneracy  of 
a  given  substate  have  been  related  to  the  quantum  numbers  i  and  N.  The 
average  degeneracy  of  substates  within  the  critical  interval  has  been 
determined  from  these  relations  to  first  order  in  1/N.  By  using  the 
energy  levels  and  average  degeneracies  computed  in  this  way,  the  ioniza¬ 
tion  and  recombination  rates  have  been  minimized  with  respect  to  the 
critical  energy  level.  The  resulting  expressions  permit  determination  of 
the  critical  state  for  any  (effective)  one-electron  atom  whose  ionization 
potential,  effective  ground  state  quantum  number,  and  maximum  orbital 
quantum  number  in  the  interval  between  ground  state  and  the  state  one 
principal  quantum  number  above  ground  are  known.  This  procedure  con¬ 
trasts  with  Byron's  method  of  determination  of  the  critical  state.  In  the 
latter,  energy  levels  and  degeneracies  of  all  states  of  interest  must  be 
used  in  a  graphical  procedure  which  must  be  repeated  for  each  separate 
electron  temperature.  The  critical  states  of  the  present  theory  agree 
within  10  percent  with  Byron's. 

Using  the  analytical  expressions  developed  for  the  critical  state  be¬ 
havior,  the  collisional  recombination  coefficients  for  hydrogen,  cesium, 
and  potassium  have  been  presented  for  electron  temperature  ranging  to 
SSOO^K.  Within  the  resolution  of  present-day  experiments,  the  basic 
theory  is  well  confirmed.  Comparison  with  other  theories  and  evaluation 
of  the  limits  of  validity  of  the  mathematical  approximation  indicates,  how¬ 
ever,  that  for  sodium,  lithium,  and  potassium  the  present  theory  under¬ 
estimates  recombination  rates  at  high  temperature  by  about  50  percent 
(well  within  the  experimental  resolution);  whereas,  for  cesium  and 
rubidium,  the  present  approximations  are  quite  adequate.  Since  the 
present  theory  represents,  to  the  author's  knowledge,  the  only  analytical 
study  of  alkali  collisional  ionization  and  recombination  to  date  (the  other 
theories  are  numerical  in  character),  it  has  the  merit  of  providing  a 
simple  description  of  the  ionization- recombination  process  for  any  atom 
whose  quantum  mechanical  structure  can  be  reduced  to  an  effective  one- 
electron  description.  In  addition  to  this  inherent  simplicity  in  calculation, 
the  present  theory  advantageously  incorporates  the  effect  of  different 
atomic  structure  of  the  elements  under  consideration  through  the  depend¬ 
ence  of  the  critical  energy  level  on  the  simple  parameter  o.  This  theory 
is  easily  generalized  through  use  of  the  methods  of  Section  VI  to  include 
radiative  effects,  provided  radiative  transitions  are  not  the  dominant  mode 
for  traversing  the  critical  interval.  For  collision- dominant  conditions,  the 
present  theory  is  accurate  within  (at  worst)  50  percent.  For  radiation- 
dominant  conditions  the  Hinnov  and  Hirschberg  theory  is  more  applicable 
than  either  the  present  theory  or  the  various  forms  of  the  Byron  theory. 
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At  present,  truly  definitive  experiments  concerning  ionization  and 
recombination  in  collision- dominated  alkali  plasmas  do  not  exist.  With¬ 
in  the  resolution  of  present-day  experiments  (factor  of  two),  the  general 
theoretical  model  is  well  confirmed.  An  additional  test  of  the  model  can 
be  obtained  by  comparing  experimentally  observed  relaxation  in  stream¬ 
ing  plasmas  subjected  to  a  constant  electric  field  with  theoretical  pre¬ 
dictions. 

A  computer  program  is  being  developed  at  AEDC  whereby  the  non¬ 
linear  electron  equations  of  motion  can  be  solved  numerically,  using  the 
ionization  and  recombination  theory  developed  in  this  report.  Such  a 
procedure  has  the  advantage  that  it  not  only  produces  relaxation  profiles 
of  electron  temperature  and  number  density  (to  be  compared  with  exper¬ 
imental  data),  but  also  yields  the  theoretical  time  variation  of  the  critical 
state,  which  can  be  compared  with  time- resolved  spectra.  The  disad¬ 
vantage  of  this  procedure  is  relative  insensitivity  of  the  results  to  the 
ionization  and  recombination  cross  sections,  complicated  by  the  lack  of 
precise  knowledge  of  possible  electron  energy  losses  within  the  plasma. 

Within  the  restrictions  imposed  by  these  ambiguities,  preliminary 
results  of  the  calculations  indicate  that  the  experimental  data  in  Ref.  35 
on  relaxation  of  a  potassium- seeded  argon  plasma  can  be  fairly  well 
corroborated  with  theory. 

More  refined  experiments  along  the  lines  of  interest  of  the  present 
study  should  ideally  include  both  electrodynamic  and  spectral  measure¬ 
ments  at  various  locations  along  a  channel  through  which  a  high-speed 
(1000  m/sec)  low  temperature  (2000'’K  or  less)  alkali-seeded  plasma  is 
flowing.  In  order  to  avoid  low  electron  density  difficulties  the  electric 
fields  imposed  should  be  at  least  10  v/cm.  Visible  radiation  should  be 
monitored  at  least  every  5  cm. 
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Fig.  1  Ratio  o1  Excited  State  Populations  in  Saha  Equilibrium 
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Fig.  2  Alkali  Energy  Level  Diagram  Showing  Effect  of  Collisionol  Transitions 
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Fig.  3  Critical  Energy  Level  in  Equilibrium  with  Free  Electrons 
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Fig.  5  Coll  isionol  Recombination  Rates  for  Hydrogenic  Elements 
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Fig.  £  Electron-Ion  Recombination  Coefficient  os  o  Function  of  Electron  Density 
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TABLE  1 

INTENSITY  OF  ALLOWED  RADIATIVE  TRANSITIONS 

ACROSS  VARIOUS  CRITICAL  INTERVALS 

INTERVAL  TRANSITION 

OSCILLATOR  STRENGTH 

7S  -  6S  -  6Pv 

0.21 

6P«,  -  6S., 

0.81 

Cnmposhc 

0.41 

7S.  ■  —  6Pi 

X  i 

0.17 

6Pi  -  6Si, 

0.39 

'*  a 

Composite 

0-26 

6D./  -  6P./ 

0.33 

4  '1 

6Pv  -  6Sv 

/a  4 

0.81 

Composil'^ 

0.52 

61).  -  6P. 

0.01 

6Pi  —  <jSi 

0.81 

3  ‘i 

Compoaite 

O.IB 

6Dv  -  6Pv 

t  'i 

0.30 

6?u  -  6Sv 

'i 

0.39 
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0.34 

:p>-  -  6s.. 

0.02 

a  i 

Sum  of  OsL'ilUtur  Stiengths 

1.73 

7P  -  6P  TPs  -  TS. 

1+ 

•a  'a 

7Si  “  6Pj/ 

0.21 

a  o 

Compos  lie 

0.46 

TPs.  -  TSi. 

1-. 

»  'J 

TS.  -  6Pi, 

0.17 

a  a 

(^onpoHiie 

0.41 

TP.  -  aD., 

1- 

50,.  -  CP. 
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Coo^posiic 

0.45 

TP.,  -  51)., 
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5P.,  -  6P.. 

0.25 

a  a 

Compiisilc 

0.23 

TP.  -  5D.. 

0,65 

a  s 

51).,  -  6P. 
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'a  a 
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0.41 

\ 

7Pv  —  6Si 

4  fz 

0.02 

8S.,  -  6P. 

0.02 

a  a 

8S.  -  6Pi. 

0.02 

a  3 

Sun  of  Oscillator  Strengths 

2.01 
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TABLE  I  (Concluded) 


6D 

-  5D  -  7Pi. 

'i  ^ 
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4r  -  SDa, 

*1 

0.32 

4F  -  5Di/ 
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'a  'a 
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n  U 
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h  n 

0.02 

8Si/  -  6Pi., 

••a  'a 

0.02 
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TABLE  II 

EFFECTIVENESS  OF  LINE  RADIATION  VERSUS  COLLISIONAL  TRANSITIONS 
ACROSS  VARIOUS  CRITICAL  INTERVALS 


INTERVAL 

SUMMED  OSCILLATOR  STRENGTHS 

"e  X 

( 

i 

cm 

7S 

-  6S 

1.73 

1.24 

X 

10“ 

7P 

-  6P 

2.01 

2.18 

X 

lo” 

6D 

-  5D 

3.29 

2.13 

X 

10“ 

as 

-  7S 

1,68 

2,88 

X 

lo” 

8P 

-  7P 

1.19 

5.18 

X 

10“ 
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